Spin Current in p-wave Superconducting Rings 
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A formula of spin currents in mesoscopic superconductors is derived from the mean-field theory 
of superconductivity. Spin flow is generated by spatial gradient of d which represents a spin state 
of spin-triplet superconductors. We discuss a possibility of circulating spin currents in isolated p- 
wave superconducting rings at the zero magnetic field. The direction of spin currents depends on 
topological numbers which characterize spatial configurations of d on the ring. 
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The generation and control of spin polarized currents 
have been the main topics in spintronicsi^. Spin devices 
have been originally proposed on ferromagnets or half- 
metals because electric currents in such materials basi- 
cally carry spins at the same time. To realize more effi- 
cient spin devices, hybrid structures of ferromagnets and 
semiconductors^ may be more advantageous. In this re- 
search direction, ferromagnetic semiconductors with high 
Curie temperatures are desired to achieve high spin in- 
jection rates into semiconductors. Ferromagnetism, how- 
ever, is not always necessary to generate spin currents. 
In a recent theoretical study, the dissipationless spin cur- 
rent due to the quantum geometric phase was proposed 
in nonmagnetic semiconductors 4 . In contrast to ferro- 
magnetic metals and semiconductors, superconductors 
are not major materials in spintronics. This might be 
because Cooper pairs do not carry spins in conventional 
s-wave superconductors. Cooper pairs in p-wave super- 
conductors, however, have the spin degree of freedom. 
Thus spin active transport can be expected in spin-triplet 
superconductors^. 

Spin states of Cooper pairs are described by three com- 
ponents of d. Throughout this paper, vectors in spin 
space are denoted by • — . Properties of d are similar to 
those of magnetic moments to in ferromagnets because 
both of them characterize spin polarization. Here we 
summarize three important differences between d and to. 
Firstly the time-reversal symmetry (TRS) is broken in 
ferromagnets, whereas it is preserved in unitary super- 
conducting states ( i.e., d x d* = 0). It is known that 
breaking down of the TRS is necessary to generate elec- 
tric current in equilibrium. Secondly to is an observable 
in ferromagnets, whereas d itself is nonobservable. This is 
because d is a part of the wave function of Cooper pairs. 
The fractional-flux states were discussed 6 in twisted crys- 
tals! based on this property. Finally m is parallel to spin 
polarization, while d points the perpendicular direction 
to spin polarization of Cooper pairs. 

In this paper, we derive a formula of spin currents 
based on the mean-field theory of superconductivity. To 
overview basic properties of spin currents, we first ap- 
ply the obtained formula to Josephson junctions of spin- 



triplet superconductors. Spin currents are represented 
by the Andreev reflection coefficients, which implies that 
Cooper pairs carry spins. We mainly discuss the circulat- 
ing spin current in isolated p-wave superconducting rings 
by using the obtained formula. The direction of spin 
current depends on spatial configurations of d which are 
characterized by topological winding numbers. Through- 
out this paper, we take the unit of K = kg = c = 1, where 
ks is the Boltzmann constant and c denotes the speed of 
light. 

Electronic states in superconductors are described by 
the Bogoliubov-de Gennes equation, 
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where D r = V r — ie A r , • ~ • indicates 2x2 matrix describ- 
ing spin space, <to is the unit matrix and /i is the Fermi 
energy. In uniform superconductors, pair potentials are 
given by 



A(r,r') =J-VA fe e ife -(— '), 



Afc =i d(k) • <7<72e z 



(3) 
(4) 



where &i with i = 1 , 2 and 3 are the Pauli matrices, and 
tp is a macroscopic phase. Using a conservation low of 
spin density, we derive a formula of the spin current, 
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where & UIn (r,r') is the Matsubara-Green function of 
Eq. in the 4x4 matrix form, and Tr is carried out 
over the Nambu x spin space. When a/2 is replaced by 
— e&o in Eq. we obtain the formula of the Josephson 
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electric current (J e )&S. Eq. J5J) is a general expression 
of spin currents which can be applied to any supercon- 
ducting systems. From Eq. (JjJJ, it is possible to derive a 
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FIG. 1: Fhe Josephson junction under consideration is in (a). 
In (b)-(d), schematic pictures of p-wave rings are illustrated. 

spin current formula in superconductor / superconductor 
junctions. When triplet superconductors are in unitary 
states, spin currents are calculated from the Green func- 
tion of junctions^, 



T 



lis » O „ - <T* 

_jAL I+ ai-+a 1 A 2#1+T 



L,- 



(6) 



where fl Li ± = y u% + \d Lj ±\ 2 , A L:± = i d L ,± ■ aa 2 , 

6l,± — d/,(±fc,p) are the vectors in the left supercon- 
ductor. In the momentum space, k and p are wave num- 
bers on the Fermi surface in the direction of currents and 
in the directions transverse to currents, respectively. The 
electric Josephson current is also given by a/2 — > —e&o in 
Eq. © 8 . Most important feature of Eq. © is that spin 
currents are represented by the two Andreev reflection 
coefficients of a quasiparticle incident from a supercon- 
ductor on the left hand side, (i.e., Si and a 2 ). This fact 
implies that Cooper pairs carry spins. 

To overview basic properties of spin currents, let us 
briefly discuss spin currents in a Josephson junction in 
the clean limit as shown in Fig. [Tfa). We assume that 
dh(R) ar e real vectors satisfying dum^ — d^rm _ = 
d>L(R) f° r simplicity. We also assume that amplitudes 
of d in two superconductors are identical to each other, 
(i.e., \dt\ = \da\ = \d\). The directions of dn is ori- 
ented by a from dh ■ In the absence of a potential barrier 
between the two superconductors, we first calculate the 
Andreev reflection coefficients by solving Bogoliubov-de 
Gennes equation. Electric and spin currents of junctions 



in Fig. [Ja) result in, 
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where rj± = (tpia)/2, n — dji x d^/\dji x dt\. Current- 
phase relations of charge are always the odd function of 
ip irrespective of a. At a = 0, the current-phase relation 
becomes J e oc sin(</?/2) at T = for — n < ip < tt as 
it is in usual short junctions in the clean limiliiS. From 
Eq. (JSJ, it is confirmed that the current-phase relation of 
spin are always the even function of ip. Spin currents are 
allowed even at <p = (in the presence of TRS). Actually 
at (p = 0, we find 



sin(a/2) 
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for — 7r < a < 7T, where T c is the critical temperature. 
It is also confirmed that spin currents vanish at a = ±7r 
and 0. Spin flows are possible when d in the two su- 
perconductors are not aligned in parallel or antiparallel. 
The calculated results in Eq. (JSJ) seem to be consistent 
with those obtained by the quasiclassical Green function 
theoryii. The spatial part of pair potentials affects am- 
plitudes of spin currents through the integration of |d| 
on the Fermi surface as shown in Eq.JSJ). The difference 
of pairing symmetries such as p- and /-wave symmetries 
is not so important in weak links. In the presence of a 
potential barrier between the two superconductors, how- 
ever, properties of spin currents may depend on the spa- 
tial part of pair potentials. Amplitudes of pair potentials 
at the interface may deviate from their bulk values in 
real materials. Thus it needs to determine amplitudes 
of pair potentials in a self-consistent way in order to dis- 
cuss electric and spin currents in realistic junctions. Such 
self-consistency of pair potentials is not considered here 
because purpose of this part is to explain basic proper- 
ties of spin current by using the analytic representation 
in Eq. JTDJl. 

The results in Eq. I|1(J|I indicate that spin currents are 
caused by the spatial gradient of d. Generally speak- 
ing, directions of d are fixed to underlying lattices of 
superconductors because they are determined by spin 
anisotropy due to the spin-orbit coupling. Thus the non- 
homogeneity of d are not expected in usual bulk samples 
and thin films. The situation, however, would be changed 
in topological crystals^ as shown in Fig. ^b)-(c). In (b), 
no spin current is expected because d is homogeneous on 
the ring. In (c), however, it is easily found that d is ro- 
tated by 2-7T while it circles the ring. The pair potential 
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in Fig. [Jc) can be described by 

A e (k,p) = ie r {6) ■ da 2 d(k lP ), (11) 

where d(k,p) is the spatial part of the pair potential, 
e r is the unit vector in the r direction. As shown in 
Fig. IT^ d). wave numbers in 8 and (r, z) directions are 



denoted by k and p, respectively. We first calculate the 
Green function of the two-dimensional ring as shown in 
Fig. Ut c ) with Eq. Ijlljl • When the radius of ring is fixed 
at r, the Green function for z = z' and 8 > 8' are given 
by 
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where 7 = /c_p/2mr 2 , i± = (cr ± <73)/2, f2^r = 
^/Id^l 2 - (iwn ± 7) 2 with v = + or — , d± = d(±k,p), 
and kp is the dimensionless Fermi wave number. On the 
way to Eq. (|12|l . we assume a relation 7 < |<i+,-| <C kpj- 
By substituting Eq. I|12|) into Eq. JSJ), spin currents of 
the ring in Fig. ^ (c) become at the zero temperature 

Js = {-no) (|) \ 7 N C: (13) 

where e z is the unit spin vector in the z direction, rig is 
the unit vector in real space in the 8 direction, and N c 
is the number of propagating channels on the Fermi sur- 
face. We have also assumed that d + = <i_ = d. Under 
the condition 7 < the final expression of spin cur- 
rent in Eq. (|13f) depends only on shapes of rings such as 
diameters and heights. Spins polarized in the +z direc- 
tion flow in clockwise. In other words, spins polarized in 
the — z direction flow in counterclockwise. In the pres- 
ence of rotation in d, circulating spin currents flow in 
equilibrium in p-wave superconducting rings because of 
the non-commutativity of spin algebra. In small ferro- 
magnetic rings, it has been pointed out that the non- 
commutativity of spin algebra is a source of the circulat- 
ing electric current^. 

Directions of spin currents depend on configurations of 
d. In Fig. 13 we show top views of several p-wave rings 
with various configurations of d. At first sight, spin cur- 
rents in the four rings in (a) and (b) seem to be different 
from one another. The microscopic calculation, however, 
show that spin currents in the four configurations are 
the same with one another. The d in the left ring in (a) 
point opposite directions to those in the right ring in (a). 
The flip of d can be described by the 7r-phase shift in 



macroscopic phase (i.e., d — » —d = de lrr ). It is evident 
that the uniform phase shift does not affect physical val- 
ues. The same argument is valid between the two rings 
in Fig- Etb) . Configurations of d in (a) and those in (b) 
can be related to each other by the uniform rotation of 
d by 7r/2. These spatial configurations of d on rings are 
characterized by a topological winding number (N). In 
the four figures in (a) and (b), d is rotated by 2ir while 
8 increasing from to 2ir, (i.e., N = 1). Spin currents 
polarized in +z direction circulate clockwise for N = 1. 
To change the direction of spin currents, we have to con- 
sider d as shown in Fig. |2Ic), where the winding number 
of d becomes N = — 1. It is easy to confirm that d suf- 
fers the rotation by — 2ir while 8 increasing from to 2ir. 
All rings in (a) and (b) show uniform configurations of d 
when we cut a ring at a certain place and deform the ring 
into the strip. In (c), however, the nonhomogeneity in d 
persists even when the ring is deformed to the strip. A 
possibility of the d in (c) is questionable in real supercon- 
ductors because the nonhomogeneity in d costs energy. 
Therefore spin-up and -down symmetry of Cooper pairs 
is violated in their directions of flow. Spins polarized +z 
(—z) direction always flow in clockwise (conterclockwise) 
on ring crystals of superconductors. This property is a 
new feature of spin-triplet superconductors. 

In summary, we have derived a formula of spin cur- 
rents in spin-triplet superconducting systems. Flow of 
spin in equilibrium is described by the Andreev reflection 
and is possible when d has spatial gradient. In p-wave 
superconducting rings, directions of spin currents are de- 
termined by winding numbers which characterize spatial 
configurations of d. The spin-up and -down symmetry in 
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FIG. 2: Various configurations of d on topological ring crys- 
tals. The direction of spin currents is the same in four rings 
in (a) and (b). In (c), spins flow opposite direction to those 
in (a) and (b). 



Cooper pair is violated with respect to their directions of 
flow; this is a new feature of spin-triplet superconductors. 
In this paper, we focus on spin-triplet unitary states in 
the clean limit. Effects of nonunitary states and those of 
midgap Andreev resonant states^ on the spin transport 
are important open questions. 
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